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The qubit entanglement induced by quasiparticle excitations in the Heisenberg spin chain and 
its relationship to the Bethe Ansatz structure of the eigenmodes is studied. A phenomenon called 
entanglement quenching, which suppresses eigenstate entanglement, is described and shown to be 
mediated by Goldstone magnons. Scattering states are characterized by short-range entanglement, 
and never exhibit entanglement at the longest range. In contrast, bound states have complex long- 
range entanglement structures. 
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I. INTRODUCTION 

The importance of quantum entanglement as a tech- 
nological resource forquantum computing and quan- 
tum communication |2J, as an experimental requirement 
for Bell's inequality violations [3] and other founda- 
tional tests of quantum mechanics, and as a central sub- 
ject in the mathematical analysis of highly correlated 
states is now widely recognized. A fundamental the- 
oretical problem is to understand the complex entan- 
glement structure naturally present in systems as di- 
verse as the Rindler vacuum the photons produced 
through spontaneous parametric downconversion , and 
the exchange-coupled spins of a magnet ■ Much of the 
work on this problem to date has focused on ground state 
entanglement and its relation to quantum phase transi- 
tions y| , or on thermal state entanglement and its mod- 
ulation by magnetic fields 0, 0- In contrast, little is 
known about entanglement in excited states, or how this 
entanglement is related to quasiparticle interactions. The 
purpose of this paper is to address these latter issues by 
examining entanglement in a model system of interacting 
qubits, the Heisenberg spin chain. 

The Heisenberg spin chain (HSC) is a one-dimensional 
lattice of spin-1/2 particles, with nearest-neighbor spins 
coupled by the exchange interaction. The HSC Hamilto- 
nian is 



Hhsc — J ^ S< • Sj + i, 



(1) 



i=i 



where Sj = (Sf , Sf , Sf) is the spin operator associated 
with the spin at the ith site of the lattice; its Carte- 
sian components obey the usual angular momentum alge- 
bra. Periodic boundary conditions are assumed, so that 
Syv+i = Si. J is the exchange coupling constant. States 
in the 2 Ar -dimensional Hilbert space Ti. of the model can 
be expanded in the "standard" basis of states of the form 
■ ■ ■ , k), where the listed spins i,j,...,k are aligned 
in the +z-direction ('up') and the remaining spins are 
antiparallel ('down'). |0) denotes the state with all spins 
down. The HSC and its anisotropic relatives have been 
the subject of considerableprevious research in quantum 
information theory 0,0,13. 

The HSC model was originally solved by a method now 



termed the Bethe Ansatz (BA) [T0| . In the physical in- 
terpretation of this solution eigenstates are comprised of 
interacting quasiparticle excitations called magnons; the 
ath magnon is characterized by a pseudomomentum k a , 
and there is a phase (j) a b between magnons a and b. We 
would like to understand how the entanglement between 
the spins of the HSC, viewed as qubits, depends on the 
pseudomomenta and phases of the magnons. Unfortu- 
nately, as we will see, there is in general no simple rela- 
tionship connecting qubit entanglement to these parame- 
ters. Nevertheless, the Bethe Ansatz formalism provides 
information that is not available from a direct numeri- 
cal diagonalization of the Hamiltonian (although the BA 
equations themselves often must be solved numerically), 
and this information has much to tell us about the struc- 
ture of qubit entanglement in the spin chain. Because the 
BA formalism is essential to the discussion, it is reviewed 
very briefly in the following section. 

The entanglement between two spins will be quantified 
by calculating their concurrence , a measure of the in- 
separability of two-qubit pure or mixed states commonly 
used in quantum information theory. Other metrics, such 
as the entanglement of formation, can be calculated di- 
rectly from the concurrence. Because the HSC Hamil- 
tonian Hhsc commutes with the total z-spin operator 
S z = J2iLi &i> the reduced density matrix (RDM) for 
any two spins in the states of interest always takes the 
form 



a 
(3 7 
7* S 
e 

For such a density matrix, the concurrence is 



C = 2max(0, | 7 | 



(2) 



(3) 



Concurrence ranges from zero, for a separable state, to 
one, for a maximally entangled state. Because the states 
of the HSC which will be examined are translationally 
invariant, the concurrence depends only on the separa- 
tion between the two spins. The concurrence between 
nearest-neighbor spins will be denoted C±, between next- 
nearest- neighbor spins C2, etc. 
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II. BETHE ANSATZ SOLUTION OF THE HSC 

The Bcthc Ansatz consists of the hypothesis that the 
coefficients of the expansions of the eigenstates of the 
HSC Hamiltonian (Eq. ^) with respect to the standard 
basis have the form of a permutation-invariant sum of 
n! exponential terms, where n is the number of spins in- 
verted with respect to the reference state |0) ^j- For ex- 
ample, the two-magnon eigenstates are assumed to have 
the form 

\ip(k 1) k 2 ,<l>i2)) = v(ki,k2,<h.2) Y (4) 

mi <m 2 

^ e ik 1 m 1 +ik 2 m 2 +icf> 12 / 2 _|_ e ik 2 mi +iki m 2 -1^12/2 j m^j 

where v is a normalization constant. Substituting this 
state into the Schrodinger eigenvalue equation with 
the HSC Hamiltonian yields the constraints (the Bethe 
Ansatz equations, or BAE) which k\, k 2 , and <pi 2 must 
satisfy for the hypothesis to be correct: 



(Eq. are properly normalized physical states (but not 
HSC eigenstates). 

In bound states the magnons have complex conjugate 
pseudomomenta: fci,2 = u ± iv. The imaginary com- 
ponent v causes the probability distribution for the sep- 
aration of the two inverted spins to be maximal when 
the spins are adjacent, and to decay exponentially as the 
separation between them increases. Due to this exponen- 
tially tight binding of the inverted spins, bound states 
exhibit entanglement behavior quite different from that 
of the scattering states, which (as will be shown) is con- 
trolled entirely by the binding parameter v. The bound 
states are of two types, both of which can usefully be put 
into forms independent of the Bethe phase. The cosh- 
type bound (CB) states can be written as 



|VcM)) = v c {v) Y e iu ^ +m ^ (8) 

mi <m 2 
N 

cosh [v( — - [m 2 - mi) )] \mi,m 2 ), 



2 cot 



'11 
2 

Nkl - 012 

Nk 2 + 0i2 



= cot 



— cot 



(5) 



2ttAi, 
2ttA 2 . 



Note that under the interchange fej «-> k 2 , (f>\ 2 changes 
sign, so that the wavefunction coefficients in Eq. 0] are 
permutation-invariant, as desired. \\ and A2 are inte- 
gral quantum numbers between and N — 1 whose al- 
lowed values are prescribed by an arcane set of rules ^3 • 
k\ and k 2 are interpreted as the pseudomomenta of two 
magnons, and 4>i 2 is a phase (the Bethe phase) between 
them. 

The BA leads to a natural classification of two-magnon 
states as either scattering states or bound states. In scat- 
tering states each magnon has a real pseudomomentum 
between and 2ir, and the Bethe phase is real and be- 
tween and 7r if one chooses k\ <k 2 . For such states the 
state vector (Eq. [3J can be written in the form 

Wfci,fca,&2)> = 2i/(Mia) Y ( 6 ) 

mi<m 2 

e ^(mi+m 2 )/2 1 (fc(m2 _ + ^ m2)j 

where k = k 2 — k\ is the relative pseudomomentum, 
K = k 2 + k\ is the total pseudomomentum, and the nor- 
malization constant v{k, <j>n) is 



v 2 {k,4>\ 2 ) 
N cos(fc - 



= N(N — 1)- 
cf> 12 )-(N- 



1) cos0 12 — cos(kN — 4>i 2 ) 



1 — cos(fc) 



(7) 



This expression remains correct in the limit k — > 0. Even 
when the BAE are not satisfied, states of the above form 



where 



uc(v) 



4 sinhf 



, (9) 



_7V[sinh((A- l)v) + (N — 1) sinhw] 
while the sinh-type bound (SB) states can be written as 



\ips(u,v)) = u s (v) Y 



iu{m\ -\-rn2) 



l/ S [V) > ( 

mi <m 2 

N 

sinh [v( — - (m 2 - mi) )] \mi,m 2 ), 



(10) 



where 



u s (v) 



4 sinh v 



N [sinh ((N - l)v) - (N - 1) sinh?;] 



(11) 



More generally, when there are n magnons, the n pseu- 
domomenta and n{n — l)/2 phases satisfy a set of n lin- 
ear and n(n — 1)/2 transcendental coupled BA equations. 
Multimagnon states can consist of all scattering magnons 
with real pseudomomenta, but more exotic states also 
appear. These include "wavecomplexes" (to use Bethe's 
terminology [Toj^. which are solitons in which a group 
of magnons are all mutually bound, and mixtures of one 
or more wavecomplexes with scattering magnons. Note 
that instead of choosing as the reference state |0), the 
state with all spins down, one could equally well have 
chosen the state with all spins up, and statements which 
hold for n-magnon states can therefore be translated into 
equivalent statements for (JV — n)-magnon states. 
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III. ENTANGLEMENT QUENCHING 
A. Two-magnon entanglement quenching 

The zero-magnon eigenstate of the HSC is the refer- 
ence state |0); it is of course a pure product state. The 
excitation of a single magnon changes this dramatically. 
The one-magnon eigenstates have the expected BA form 



0.2 



l^(fci)) 



N 

E 

rn — 1 



e ikim \m) 



(12) 



where fci = 2ttXi/N and Ai is a quantum number that 
takes integral values between and N — 1. Remarkably, 
the spins of the HSC are now equientangled: the concur- 
rence between any two spins is 2/N irrespective of the 
pseudomomcntum of the magnon or of the separation of 
the spins on the chain. 

What is the effect on entanglement of the excitation 
of a second magnon? To begin, we consider the spe- 
cial case where the second magnon has zero pseudo- 
momentum. Because of the magnon dispersion relation 
Ek ~ 1 — cos(fc), no energy is required to excite such a 
magnon: it is "massless" . Such modes can be considered 
as arising from the spontaneous breaking of the global 
SU(2) symmetry of the HSC Hamiltonian (Eq. by 
the choice of the reference state 10), and hence will be 



referred to as Goldstone magnons 

From the two-magnon BA equations (Eq. |SJ, we see 
that when one of the pseudomomenta vanishes, the Bethe 
phase 0i2 must also vanish. The HSC is now in the state 

IV (Ai =O,fc 2 ,0 12 =0)> = 2i/(A a> 0) E 

mi<m 2 

e ik 2 (m 1+m2 )/2 CQS ( fc2 ( m2 _ mi )/2) \mx,m 2 ). (13) 

We can visualize the entanglement in this state by plot- 
ting the concurrence between two spins as a function 
of the relative pseudomomentum k = fc 2 , treating k as 
a continuous parameter, as in Fig. ^ These states are 
eigenstates precisely when k = 27rA 2 / 'N (A 2 = 0, . . . , N — 
1), which is also the condition for translational invari- 
ance, and so only the separation of the spins matters at 
the points of interest. When k = 0, all qubit concurrences 
are equal with Q w 0.21. Strikingly, all other eigenstates 
correspond to exact zeros of the concurrences Q. The ex- 
citation of a Goldstone magnon has completely quenched 
the qubit equientanglement, except in the special case in 
which the original magnon had zero pseudomomentum 
as well. 

It is possible to prove analytically that this must al- 
ways occur whenever a Goldstone and a non-Goldstonc 
magnon mix. Choose two arbitrary spins p < q of the 
HSC and trace out the rest: this yields the reduced den- 
sity matrix p{p,q), whose elements are (in the notation 
of Eq. EJ: 



C1 



. 1 




0.2 



C 2 



. 1 



0.2 



c 3 



. 1 



7t 

3 
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3 
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FIG. 1: Concurrence between spins 1 and 2 {top), 1 and 3 
[middle), and 1 and 4 {bottom) in a two-magnon scattering 
state \ip{0, k,Q)) of the six-spin Heisenberg chain. The rela- 
tive pseudomomentum k is treated as a continuous variable 
parameterizing states of the scattering form (Eq. ; 0i2 is 
set to zero. Eigenstates occur at k = 2n\2/ 'N , A2 = 0, . . . , 5, 
and are indicated by ticks on the fc-axis. 



a(p, q) 
0(p,q) 

lip, q) 
t(p, q) 

where 



N(N - 


-2) 


$(p, q) 


T 


2{N- 


4) 


N{N - 


-2) 


1 




N(N - 


-2) 



c 2 (A 2 



(14) 

(A-2-2c 2 (A 2 ))(15) 



N(N — 2) 
exp(TTi(p-q)\ 2 /N)c{\ 2 ) (16) 



((A-2)(A-4) + 4c 2 (A 2 ))(a7) 



c(A 2 ) = cos 



?rA 2 (p- q) 
N 



(18) 



and A 2 € {1, . . . , N — 1} (these matrix elements are not 
correct when A 2 = 0) . The concurrence between spins p 
and q, from Eq. |3l is 



C(p g)- A ^ x 
C[P ' q) - N{N-2) X 

max (0, N - 4 - y/(N - 2){N - 4) + 4c 2 (A 2 )) (19) 



But 



^(iV-2)(iV-4)+4c 2 (A 2 ) > y/(N-2)(N-4) 

> A- 4 (20) 
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for N > 4 and hence C(p, q) is identically zero. 

Because the two-magnon state under consideration 
(Eq. 113(1 is pure, the entanglement of formation Sf for 
the reduced density matrix p(p, q, r, . . .) of any subset 
p,q,r, . . . of the spins of the HSC coincides with the bi- 
nary von Neumann entropy of that density matrix: Sf — 
-Tr[p(p,q,r,...)lgp(p, q,r, ...)], where \g(z) = log 2 (x). 
In particular, the entanglement of formation for the state 
of a single spin p is simply 

*<~®*G)-H)*H)-™ 

This is nonzero for any finite N > 2. These results yield 
the somewhat paradoxical conclusion that in states con- 
taining both a Goldstone and a non-Goldstone magnon, 
any specific spin is not entangled with any other spin in- 
dividually, although it is entangled with all other spins 
collectively. 



B. Multimagnon entanglement quenching 

Suppose that a third magnon is excited on a HSC 
initially in a two-magnon entanglement-quenched state 
(Eq. 113(1 . Will this create qubit entanglement? Again we 
can begin with the special case in which the additional 
magnon is also a Goldstone excitation. Three-magnon 
state vectors are somewhat difficult to manipulate ana- 
lytically. They can however be generated with arbitrary 
precision using a numerical implementation of the Bethe 
Ansatz, and the partial trace required to find the two- 
spin RDM and the concurrence can then be done numer- 
ically with no loss of accuracy. All possible concurrences 
in all states describing the simultaneous excitation of two 
Goldstone and one non-Goldstone magnons for Heisen- 
berg spin chains of lengths from N = 6 to N = 50 were 
calculated in this manner and found to be zero. Thus 
a second Goldstone magnon cannot revive quenched en- 
tanglement. 

A different situation arises when two non-Goldstone 
magnons are present. Such states have complex entan- 
glement structures, as will be discussed later. Table Q] 
shows the effects of the excitation of a Goldstone mode 
on qubit entanglement in the two-magnon scattering and 
bound states of an N = 6 ring (in this case the BAE can 
be handled analytically). In five of the states entangle- 
ment is completely quenched. In two of the bound states 
(quantum numbers Ai = A 2 = 1 or 5) the qubit entan- 
glement is reduced, but not eliminated, while in two of 
the scatting states (Ai = 1, A 2 = 5 and Ai = 2, A 2 = 4), 
entanglement is actually generated at the longest range 
(i.e. C3). Nevertheless, the total entanglement, as mea- 
sured by the sum of the concurrences Q, has decreased. 
This phenomenon is confirmed by numerical studies of 
other rings: excitation of a Goldstone magnon will occa- 
sionally generate entanglement between a formerly unen- 



tangled pair, but the total entanglement of the HSC is 
always reduced. 

Finally, we consider entanglement quenching in the 
four-magnon case. Suppose a second Goldstone magnon 
is excited on a ring with two non-Goldstone and one 
Goldstone magnons. Numerical examination of all pairs 
of spins for Heisenberg spin chains of lengths from N = 8 
to N = 20 show that all concurrences vanish in such 
states. Concurrence is also absent from these chain in 
all states with three Goldstone and one non-Goldstone 
magnons. These results, and those of the preceeding 
section, support the following hypothesis: the addition 
of a Goldstone excitation to a state containing one or 
more non-Goldstone magnons always reduces the total 
qubit entanglement in the HSC; and when the number 
of Goldstone magnons equals or exceeds the number of 
non-Goldstone excitations, there is no qubit entangle- 
ment whatsoever. 



IV. THE PURE GOLDSTONE SECTOR 

As remarked previously and illustrated in Table. ^ the 
excitation of second Goldstone magnon from a single- 
Goldstone state reduces but does not completely quench 
qubit entanglement. Thus from the standpoint of en- 
tanglement behavior, the pure Goldstone sector is fun- 
damentally different from the mixed sector. As further 
Goldstone magnons are excited, the qubit concurrence 
continues to decrease without vanishing. There is no 
length scale associated with these Goldstone magnons 
(because their pseudomomentum is zero), and the HSC 
qubits remain cquicntanglcd in all states comprised en- 
tirely of Goldstone excitations. This situation can be 
studied analytically. The general formula for the concur- 
rence between any two qubits in an n-Goldstone state of 
an A-qubit ring is 
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TABLE I: Partial quenching of entanglement. Ai and A2 
are the quantum numbers of two non-Goldstone magnons on 
a six-spin HSC. Concurrences d are shown without [cen- 
ter column) and with {right column) an additional Goldstone 
magnon. The first six entries are scattering states, while the 
next three are bound states. 
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FIG. 2: Concurrence between any two qubits in a iV-spin 
Heisenberg chain as a function of the number n of Goldstone 
magnons present. 



C(n,N) = 



N(N- 1) 
[n(N — n) — y/ n(n 



l)(JV-n)(JV-n-l)].(22) 



This result has been derived elsewhere, and its implica- 
tions for thermal entanglement have been discussed [l3j . 
When n = or n = N, the eigenstate corresponds to 
all spins pointing down (the state |0)) or all spins point- 
ing up, and hence is factorable. Maximal entanglement 
occurs when a single magnon is excited with respect to 
either of these states, i.e. n = lorn = JV-l. Exci- 
tation of a second Goldstone magnon produces a large 
relative decline in entanglement, but thereafter qubit en- 
tanglement is only weakly dependent on the number of 
Goldstone magnons excited, as shown in Fig. [21 



FIG. 3: (right) Contour plots of concurrences between spins 1 
and 2 (Ci), spins 1 and 3 (C2), etc., for two-magnon scattering 
states (Eq. |SJ of an 8-spin HSC, as functions of the relative 
pseudomomentum k and the Bethe phase <f> = (j>i2- White 
regions denote zero concurrence; contours are evenly spaced 
at intervals of 0.05 (top two graphs), 0.03 (third graph), and 
0.01 (bottom graph). Values of (k, <f>) corresponding to energy 
eigenstates are indicated by dots. Small dots denote a single 
eigenstate. The six large dots denote two eigenstates with the 
same (k, <f>) values and hence the same entanglement; such a 
pair need not be energetically degenerate. (Note that near 
k » 1, 4> « 2.5 a single state and a double state nearly co- 
incide.) The point at (0,0) is the eigenstate containing two 
Goldstone magnons. The seven eigenstates with one Gold- 
stone magnon form a row on the c6-axis. In the C4 plot the 
dotted lines (for example, bisecting the central "islands") are 
zero-entanglement boundaries between two entangled regions. 
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FIG. 4: Statistical characterization of entanglement in two-magnon states for the N = 40 HSC. Left column: scattering 
eigenstates. Center column: cosh-type bound eigenstates. Right column: sinh-type bound eigenstates. The upper graphs 
show the maximum concurrence at a qubit separation n among all eigenstates in the indicated population, and the median 
concurrence of the subpopulation of states with nonzero C n . The middle graphs show the percentage of two-magnon eigenstates 
with nonzero C„ as a function of n. The lower graphs show the percentage of two-magnon scattering states states with exactly 
v values of C n not equal to zero. 



V. THE GOLDSTONE-FREE SECTOR 

A. Two-magnon scattering states 

The previous sections have characterized qubit en- 
tanglement in the pure Goldstone and the mixed 
Goldstone/non-Goldstone sectors of the HSC. We now 
consider the pure non- Goldstone sector, beginning with 
the two-magnon states. For scattering states (Eq. |SJ) it 
is easy to show that the concurrence depends only on 
the relative pseudomomentum k and the Bethe phase t/>, 
and not on the total pseudomomentum K . We can plot 
the concurrence as a function of these variables uncon- 
strained by the BAE (Eq. ISJ. This is shown in Fig. |3 for 
the N — 8 case, with the eigenstates (those parameter 
values which do satisfy the BAE) indicated by dots. The 
fifteen zero-Goldstone states form an irregular wedge in 
the upper half-plane. The double-Goldstone state and 
the single-Goldstone states appear as a row on the cf>- 



axis. As in Fig. ^ the positioning of these latter states 
is striking: although always unentangled, they often lie 
near or on the boundaries of regions of entanglement, so 
that a slight perturbation of the Bethe parameters would 
yield an entangled state. There is no relation between 
an eigenstate's entanglement and its eigenenergy (which, 
unlike its concurrences, depends on the total pseudomo- 
mentum) , and, as the complexity of the concurrence di- 
agrams suggests, there is no simple relationship between 
entanglement and the Bethe parameters. 

Although qubit entanglement cannot easily be calcu- 
lated directly in terms of the Bethe parameters for non- 
Goldstone multimagnon states, the Bethe Ansatz classi- 
fication partitions eigenstates into populations with sta- 
tistically distinct patterns of qubit entanglement. In the 
two-magnon case, the scattering states are characterized 
qualitatively by short-range entanglement, which distin- 
guishes them from the bound states. There are a number 
of ways to make this observation quantitative. One way 
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is to ask about the maximum concurrence attained by 
any scattering state as a function of qubit separation. 
As shown in Fig. 01 for the N = 40 case, this maximum 
decreases roughy linearly (but nonmonotonically) with 
increasing qubit separation, and in fact is exactly zero at 
the longest range (C20), that is, between opposite spins 
on the Heisenberg chain. This is not simply due to out- 
lying extreme values. If one culls the total population of 
two-magnon scattering states, and calculates the median 
concurrence as a function of qubit separation n among 
those states with nonzero C n , the median of each subpop- 
ulation is typically between one-half and one-third of the 
maximum. Another measure of the short-range nature of 
entanglement is the decline in the number of eigenstates 
with nonzero concurrence as a function of distance. Over 
60% of scattering states have nearest-neighbor entangle- 
ment; this decays roughly linearly with increasing qubit 
separation, and as already noted becomes zero at the 
longest range. Most two-magnon scattering states have 
entanglement at only a few different distances, usually 
only 5 or 6; no scattering state has more than 11 nonzero 
values of C n . Similar observations hold for other values 
of AT. 

Examination of C4 entanglement in Fig. [3] (where 
N = 8) shows that all the scattering states lie in re- 
gions of zero concurrence. Similarly, inspection of Ta- 
ble shows that C3 vanishes for scattering states when 
N — 6. Analytic calculations showed that C2 was always 
zero for scattering states when N = 4 or 5, and C20 was 
seen above to vanish when N = 40. These observations 
suggested the hypothesis that concurrence at the longest 
range, CW/21, is always zero for two-magnon scattering 
states. All such states for TV = 4 up to N = 50 were 
generated numerically, and C[jv/2j always equalled zero 
as conjectured. This is yet another manifestation of the 
short-range nature of scattering state entanglement. 



B. Two-magnon bound states 

In contrast to the scattering states, the two-magnon 
bound states are characterized by long-range entangle- 
ment. The behavior of the cosh-type (Eq. [SJ) and sinh- 
type fEq. HOf) states is similar but distinct. Analogous to 
the independence of scattering state entanglement from 
the total pseudomomentum K, bound state entangle- 
ment is easily shown to be independent of the phase pa- 
rameter u; qubit concurrence is controlled entirely by the 
binding parameter v. This dependence can be plotted by 
treating v as a continuous parameter unconstrained by 
the BAE (Eq. 0, as shown in Fig. [5] for the CB states 
and in Fig. [21 for the SB states for a 10-spin HSC. In the 
limit v — ► 0, the CB state becomes the two-Goldstonc 
state and all concurrences are equal. As v increases, the 
short-range components of entanglement, C\ and C2, fall 
off, while the longer-range components increase. In this 
low-w 'anomalous' parameter regime C n increases with 
n, that is, spins actually become more highly entangled 




FIG. 5: Nearest-neighbor concurrence C\, next-nearest- 
neighbor concurrence C2, etc., for cosh- type bound (CB) 
states in the 10-spin Heisenberg chain. Two eigenstates occur 
at v = 0.258 and two at v = 1.174. An 'anomalous' regime, in 
which entanglement increases with qubit separation, extends 
to about v « 0.7. 



with increasing separation. As v continues to increase, 
C\ becomes identically zero, and never revives. C2 passes 
through a minimum and rises to become the dominant 
component, and the order of the strengths of the longer- 
range concurrences undergoes an inversion, so that en- 
tanglement now falls off with increasing spin separation. 
Asymptotically, as v — > 00, all concurrences vanish ex- 
cept for C2, whose limit will be calculated below. Simi- 
lar features occur for longer- N spin chains; in particular, 
there is always an 'anomalous' region for small v. 

The physical implications of this for eigenstates are 
determined by the distribution of solutions to the BAE 
(Eq|SJ), which for the CB states can be rewritten as 



^ 'JVu\ sinhw 

~) ~ coshi;-cos(^)' 



(23) 



where A is a quantum number. Numerically, it is found 
that the roots (values of v corresponding to eigenstates) 
of this equation always lie between iV -3 / 2 (a strict lower 
bound) and In (AT) (a strict upper bound), with one 
exception which will be described below; further, the 
roots are strongly clumped towards the lower bound. 
The result is that most allowed u-values actually lie 
in the anomalous parameter region where entanglement 
strength increases with qubit separation; this becomes 
more pronounced as N increases. This is reflected in the 
statistics of the CB state population, as shown in Fig. 0] 
In general, CB states are much more entangled than scat- 
tering states; for example, nearly a quarter of the CB 
states have nonzero concurrence at every length scale. 
The magnitude of the concurrences, however, is typically 
smaller for bound states than for scattering states. 

When > 2 is even, but not when N is odd, there is 
one solution of the two-magnon BAE in which the pseu- 
domomenta have infinite imaginary parts This is 
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FIG. 6: Next-nearest-neighbor concurrence C2, etc., for sinh- 
type bound (SB) states in the 10-spin Heisenberg chain. Two 
eigenstates occurs at v = 0.521; the singular state (Eq. 12 It 
lies at v = 00. Note the absence of both nearest-neighbor 
(Ci) and longest-range (C5) entanglement. 



the exception mentioned above, corresponding to v — > 00. 
This eigenstate with singular Bethe parameters is: 



1 



N 



l^) = 7 =^(-l) m+1 Km + l). (24) 

V ^ m =l 

It is the maximally bound state; the two inverted spins 
are always adjacent. Formally, it is a CB state when TV/2 
is even and an SB otherwise. When all spins except p and 
q are traced out, the elements of the RDM are 



a(p, q) 

i(p, q) 
e(p, q) 



TV 



p+1,9 



(25) 



S(p,q) = -(2-S p+1 , q ) (26) 



-& 



TV 

N-4 + 5 p+ i >q 



N 



(27) 
(28) 



where Sij is the Kronecker delta. The concurrence is 
therefore 



2 + 2Sn,4 
TV 



(29) 



Thus there is an unusual entanglement pattern, with 
no entanglement except between next-nearest-neighbors, 
whose concurrence is simply 2 /TV if V > 4. The four-spin 
ring is exceptional; opposite corner qubits are perfectly 
entangled. Intuitively, diagonal qubits are next-nearest- 
neighbors on both sides (in both directions around the 
ring) . It is striking that the v — > 00 limit is realized by a 
physical eigenstate. 



There are several notable differences that distinguish 
the CB and SB bound states. First, in SB states nearest- 
neighbor entanglement is always identically zero. This 
can be proven analytically, by a method resembling that 
used for the proof of Goldstone-mediated entanglement 
quenching given above, but computationally much more 
involved. Second, when TV is even, longest-range entan- 
glement Cfq/2 is also identically zero. The result follows 
immediately from calculating that the coherence 7 (no- 
tation of Eq. |2J vanishes for opposite spins, irrespective 
of the value of v. Thus SB state entanglement exhibits a 
sensitivity to the parity of V that CB states do not. Most 
significantly, however, the anomalous parameter regime 
does not in general exist for SB states; rather, at v = 0, 
the maximum concurrence typically lies at some interme- 
diate range, so that Cjv/4 (or thereabouts) is the domi- 
nant form of entanglement. As v increases, the concur- 
rence curves cross at irregular intervals. The values of v 
corresponding to eigenstates are constrained by the BAE: 



, i' Nv \ 
tanh — — = 



sinh v 



coshw — cos (^jr) 



(30) 



As with Eq. 123 the roots cluster strongly at small v val- 
ues, and the eigenstate entanglement statistics (Fig. 0J 
right column) reflect the dominance of intermediate- 
range entanglement. Thus the various classes of two- 
magnon eigenstates as determined by the Bethe Ansatz 
exhibit distinct patterns of entanglement behavior. 



C. Multimagnon states 

The entanglement characteristics of two-magnon 
bound and scattering states generalize to multimagnon 
states. For example, the quantum numbers A = 
(Ai,A2,A 3 ) = (1,3,5) specify a pure scattering state of 
three magnons with real pseudomomenta for any chain 
length TV > 6, while A = (1, 1, 1) always corresponds to a 
pure bound state, i.e. a three-magnon wavecomplex. The 
entanglement behaviors of these eigenstates are shown in 
Fig. [7| The scattering state has only nearest-neighbor 
entanglement up to V = 12; next-nearest- neighbor en- 
tanglement does not become dominant until V = 17, and 
C3 entanglement arises only at TV = 20. Thus, qualita- 
tively, the range of nonzero entanglement is always much 
smaller than the longest range possible. This behavior 
is typical of all pure scattering multimagnon states stud- 
ied so far. In contrast, the wavecomplex state initially 
(TV = 6) has only longest range entanglement (written 
as C_i, as explained in the caption for Fig. UJ; next- 
longest range entanglement C_2 arises at TV — 8, and by 
TV = 16 entanglement is present at every range. How- 
ever, the magnitude of the concurrence always increases 
as qubit separation increases, just as for a two-magnon 
CB state in the 'anomalous' regime. Thus the entan- 
glement patterns observed in the two-magnon case also 
occur in multimagnon bound and scattering states. 
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FIG. 7: Qubit concurrences in the three-magnon pure scatter- 
ing eigenstate with quantum numbers A = (1, 3, 5), and in the 
three-magnon pure bound eigenstate with quantum numbers 
A = (1, 1, 1), as functions of the HSC length N. Curves are 
drawn continuously for visual clarity, but only values at inte- 
gral N are meaningful. The notation C—i is used to denote the 
concurrence Ci jv/21+1— i, that is, C_i is the concurrence at the 
longest range (between opposite qubits on the ring), C-2 is 
the concurrence at the next-longest range, etc. The (1,1,1)- 
magnon has shorter-range concurrences (lower C_, compo- 
nents) which are not shown. 



Of particular interest is the behavior of the antiferro- 
magnetic ground state (AGS), which for even-iV chains 
is nondegenerate. In the BA formalism the AGS is spec- 
ified by the quantum numbers A = (1, 3, . . . , N — 1) 
and is a pure scattering state. It has been proven by 
O'Connor and Wootters [l^ that of all even-iV transla- 
tionally invariant states with S z — 0, the AGS has maxi- 
mal nearest-neighbor concurrence. To study this further, 
the AGS was calculated analytically for N = 4, 6 and 
numerically for N = 8, 10, 12 (work is underway to study 
longer HSCs). Two new results were obtained. First, 
all concurrences except C\ vanished; that is, the AGS 
appears (based on these few cases) to have only nearest- 
neighbor entanglement, as one might anticipate for a pure 
scattering state with the maximum number of magnons 
possible. More generally, no examples have been found of 
a chain with |_iV/2j scattering magnons having any qubit 
concurrence beyond nearest-neighbor, although this work 
is preliminary. Second, while the AGS does indeed have 



the largest nearest-neighbor entanglement, other trans- 
lationally invariant S z — eigenstates can have much 
stronger entanglement at longer ranges. For example, 
the 6-spin AGS A = (1,3,5) has C x = 0.43456, but the 
bound state A = (1, 1, 1) has C 3 = 0.70313. 

VI. DISCUSSION 

The Bcthc Ansatz provides a method for describ- 
ing eigenstates of the Heisenberg spin chain in terms 
of their pseudoparticle content. This paper has shown 
that the natural classification of eigenstates originating 
in this description predicts qualitatively and sometimes 
quantitatively the qubit entanglement of these states. 
Magnons with zero pseudomomentum, here termed Gold- 
stone magnons, effectively reduce or suppress entangle- 
ment relative to the corresponding Goldstone-free state. 
From an applied point of view this is a form of noise which 
may affect the processing, storage or coherent transmis- 
sion of quantum information in exchange-coupled qubits 
|l6j | . Pure scattering states are composed of one or more 
magnons with real pscudomomcnta. While a single such 
magnon leads to qubit equientanglement, the presence of 
additional scattering magnons favors short-range entan- 
glement: statistically, qubit entanglement becomes both 
weaker and less common in the population of scatter- 
ing states as interqubit separation increases. In particu- 
lar, the antiferromagnetic ground state seems to exhibit 
only nearest-neighbor entanglement. In contrast, bound 
states consist of two or more magnons with complex pseu- 
domomcnta. For a pair of bound magnons qubit entan- 
glement is controlled solely by the binding parameter v 
and the interqubit separation. In the large-w limit next- 
nearest- neighbor entanglement approaches 2 /TV, while all 
other forms of entanglement approach zero. This v — > 00 
limit is actually realized in even-_/V spin chains by the 
singular state. Most nonsingular two-magnon bound 
states have small values of v; in this binding parame- 
ter regime the order of the concurrence curves is inverted 
for cosh-type states, so that entanglement increases with 
qubit separation, while in sinh-type bound states entan- 
glement at intermediate lengths is favored statistically. 
Behavior similar to that of the cosh-type states is seen 
in multimagnon bound states such as the three-magnon 
A = (1, 1, 1) state. Multi-wavecomplex states, comprised 
of magnons bound into groups which are not bound to 
one another (e.g. a four-magnon state of two bound 
pairs), remain to be studied, as do mixtures of scattering 
magnons and wavecomplexes. 

The outstanding question at this point is whether this 
approach to understanding entanglement via the Bethe 
Ansatz can be extended to other models. The BA can 
also be used to solve other spin- 1/2 chains, such as the 
anisotropic deformations of the HSC (e.g. the XXZ and 
XYZ models), and it will be interesting to see if the same 
features, such as short-range entanglement in scattering 
states, reappear. More challenging is the extension of 
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this approach to Hamiltonians which include hopping, 
such as the one-dimensional Hubbard model, which are 
solved by the so-called nested Bethe Ansatz [ij} • If simi- 
lar features do reappear, as seems plausible, it may reflect 



a kind of universality arising from the BA structure of the 
solutions, allowing us to predict aspects of entanglement 
behavior across a broad range of physical systems. 
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